Macroscopic tunnel splittings in superconducting phase qubits 
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Prototype Josephsori-junction based qubit coherence times are too short for quantum computing. 
Recent experiments probing superconducting phase qubits have revealed previously unseen fine 
splittings in the transition energy spectra. These splittings have been attributed to new microscopic 
degrees of freedom (microresonators) , a previously unknown source of decoherence. We show that 
macroscopic resonant tunneling in the extremely asymmetric double well potential of the phase 
qubit can have observational consequences that are strikingly similar to the observed data. 
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Recent experiments by Simmonds et al. [l( and Cooper 
et al. |3] reveal previously unseen fine splittings in 
the transition energy spectra of superconducting phase 
qubits. These splittings are interpreted as resulting from 
coupling between the circuit's collective dynamical vari- 
able (the superconducting phase describing the coherent 
motion of a macroscopic number of Cooper pairs) and mi- 
croscopic two-level resonators, hereafter called microres- 
onators, within Josephson tunnel junctions. Microres- 
onators may be an important decoherence mechanism 
for many different superconducting qubit devices 
with broader implications for Josephson junction 
physics generally. Key questions remain however. Are all 
of the observed splittings truly a microscopic property of 
junctions? Could they instead be a macroscopic property 
of the particular circuit, or a combination of microscopic 
and macroscopic phenomena? 

In fact, macroscopic resonant tunneling (MRT) can 
produce spectral splittings in multiwell systems by lifting 
degeneracies between the states of different wells. These 
effects have been probed by Rouse et al, Friedman et 
al, and others in superconducting circuits involving 
asymmetric double wells with a few left well states, and < 
10 right-well states. MRT effects have also been demon- 
strated by Crankshaw et al. || in three-junction flux 
qubits, another system in which spurious splittings have 
been reported [9j. What is not obvious is that MRT ef- 
fects can be important for systems with extremely asym- 
metric double well potentials, like the rf SQUID phase 
qubit > that have hundreds or thousands of right well 
states. In this Letter, we analyze the phase qubit in this 
limit and show that MRT produces surprisingly complex 
observational consequences that are strikingly similar to 
some of the observed data 0,0- MRT is therefore a pos- 
sible mechanism for fine splittings in a phase qubit and 
requires further examination. 

Figure 1(a) shows the circuit schematic for an rf 
SQUID. The device is a superconducting loop of induc- 
tance L interrupted by a single Josephson junction with 
capacitance C and critical current I Cl inductively coupled 
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FIG. 1: (a) Circuit diagram for an rf SQUID phase qubit. (b) 
The device can be tuned via an inductively coupled bias line 
to give an extremely asymmetric double-well. 



to a flux-bias line. The circuit Hamiltonian is 

H = 4E cP 2 /h 2 + Ej ( 7 2 /2/3 - cos 7 - J 7 ) , (1) 

where 7 is the gauge invariant phase difference across 
the junction, p = hQ/2e is the momentum conjugate to 
7 (Q is the charge on the plates of the capacitor), (3 = 
2irI c L/$o is the modulation parameter ($0 = h/2eis the 
flux quantum) , and J = I/I c is the dimensionless current 
that is induced in the loop by the applied flux bias. The 
charging energy Ec — e 2 /2C and Josephson energy Ej = 
I c ^>q/2tt determine the regime of superconducting qubit 
behavior; for a phase qubit Ej ^> Ec- 

The shape of the circuit's potential energy function 
£/(7) depends on j3 and the bias J. For (3 < 37r/2, it is 
possible to bias the circuit so that the potential has the 
highly asymmetric double- well shape shown in Fig. 1(b), 
tuned to give a shallow upper left well with just a few 
left-localized states, denoted by \n) L , and a deep right 
well with many right-localized states, denoted by \m) R . 
Simmonds et al. |l| -motivated by a number of attrac- 
tive features including reduced quasiparticle generation, 
tunable anharmoniticity of the left well potential, induc- 
tive isolation from and reduced sensitivity to bias noise, 
and nice read-out properties-have proposed using the rf 
SQUID with an extremely asymmetric double well po- 
tential as a phase qubit |4j . 

Making a cubic approximation to the left well, we de- 
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FIG. 2: (a) Numerically computed spectrum of phase qubit when I c — 8.531 /xA, C = 1.2 pF, and L — 168 pH (/3 = 4.355). 
Energies are plotted in units of frequency. The inset shows the avoided crossing due to resonant tunnel coupling between the 
left well state \1) L and a highly excited right well state, (b) The circuit parameters give an asymmetric double well like that 
shown, (c) Wavefunctions of the k = 641 eigenstate for bias values near the avoided crossing shown in the inset, (d) Solid 
points are numerically computed sizes and locations of the splittings. Solid lines are splitting sizes derived from WKB theory. 



rive the plasma frequency for small oscillations 

WL =a; (l-/3- 2 ) 1/8 [2(J*-J)] 1/4 , (2) 
where uj — y/8E c Ej/h 2 , and 

J* = (1 - /T 2 ) 1/2 + /T 1 arccos (-/T 1 ) > 1 (3) 

is the critical bias for which the left well vanishes. Note 
that the effective critical current is I* = I C J* > I c . The 
approximate number of left-well states is 

hw L 3 V E c x 

(4) 

where AUl is the barrier height. The level spacing in 
the right well is approximately fku Rl where lur is the right 
well plasma frequency, and the number of right well states 
is approximately Nr ~ AUr/Hlur, where AUr is the 
depth of the right well. 

Figure 2(a) shows the energy spectrum as J is var- 
ied for < N L < 6 and C = 1.2 pF, L = 168 pH, 
and I c = 8.531 fiA, giving (3 = 4.355, I* = 11.659 /xA, 
and wl/27t ~ 10 GHz. These are the circuit parame- 
ters from assuming that the critical current quoted 
there is /* . To obtain the energy spectrum we diagonalizc 
the Hamiltonian in Eq. Q using a discrete Fourier grid 
representation 10], thereby obtaining a numerical solu- 
tion for the eigenvalues E^ (J) and eigenstates \k (J)) of 
the full double- well system versus the bias J. A har- 
monic approximation to the right well yields approxi- 
mately 500 states below the left well; the full calculation 
yields Nr ~ 600— 700 states, depending on the bias 

In Fig. 2(a) we define the zero of energy to be at the 
bottom of the left well. We note two different types of 
energy levels: horizontal (H) branches and near verti- 
cal (V) branches. From our definition of zero energy, 
eigenvalues corresponding to states mainly localized in 



the right well [region I of Figs. 2(a) and (b)] fall with 
increasing J, and are thus nearly vertical. The energy 
levels in region III correspond to delocalized states fully 
above the left well. The dashed line in Fig. 2(a) dividing 
regions II and III indicates the energy at the top of the 
left-well barrier. In region II, eigenstates whose energies 
lie along H branches are primarily localized in the left 
well (H ~ L). The number of left- well states at bias J 
is consistent with Nl from Eq. J2J. Eigenstates whose 
energies lie along V branches are primarily localized in 
the right well (V ~ R). Their energies fall at essentially 
the rate of the falling right well. Note that in Fig. 2(a) 
the density of right-well states is comparable to that of 
the left- well, despite N R ^> N L . 

Every apparent intersection of an H and V energy 
level in Fig. 2(a) is an avoided crossing (see inset). De- 
generacies are lifted by resonant tunneling of left-well 
states \n) L and right-well states \m) R . Left of an avoided 
crossing between k and k + 1 eigenstates we find that 
\k) = \n) L and \k + 1) = \m) R . Right of the crossing the 
states swap, becoming \k) = \fn) R and |fc + l) ~ \n) L . 
At the avoided crossing \k) = {\n) L + \m) R ) /v2 and 
\k + 1) = (\n) L — \m) R ) /y/2. Figure 2(c) shows the wave- 
functions for the k = 641 eigenstate before, at, and after 
the splitting shown in the inset in Fig. 2 (a). The distri- 
bution of splitting magnitudes along the first five energy 
branches are plotted in Fig. 2(d) as solid points. Gaps 
larger than 1 MHz arc within the resolution of recent 
experiments. Along each left-well energy branch the tun- 
nel splittings are regularly spaced with magnitudes that 
decrease exponentially with N^. We have numerically 
computed spectra for a variety of circuit parameters, in- 
cluding I c = 2 fiA and C = 0.5 pF which are comparable 
to those reported in [2j . In each case the spectrum looks 
qualitatively similar to Fig. 2 (a). We note that the pre- 
dictedgap sizes are strikingly similar to those reported 
in HIJ (~ 1-100 MHz). 
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The complex collection of energy splittings has both 
direct and indirect effects that should be taken into ac- 
count when analyzing the experimental data. Consider 
a double frequency microwave spectroscopic method, like 
that used in [1| . Microwaves of frequency wqi are applied 
to drive the — > 1 transition. Excitation of the \1) L 
state is detected with a measurement microwave pulse of 
frequency W13, which drives the 1 — > 3 transition. The 
|3}^ state's exponentially greater amplitude to be found 
in the right well compared to the |0) L and |1) L states 
allows an adjacent detection SQUID to easily detect the 
change in the qubit's flux. This method directly probes 
splittings along many of the energy branches shown in 
Fig. 2(d). Cooper et al. have introduced a new spectro- 
scopic technique that can probe deeper left wells where 
Nl > 4. This method applies a few-nanosecond current 
pulse changing the bias so that Nl > 2 by briefly tilting 
the potential adiabatically with respect to the left well 
period T L = 2tt/uj l ~ 100 ps 0. Since the measurement 
pulse moves left-well states to the right along horizontal 
(H) energy branches [see Fig. 2(a)], read-out should be 
influenced by the exponentially larger splittings present 
for smaller Nl- For example, the measurement pulse 
may move a deep well state to one of the large splitting 
degeneracies near Nl ~ 2, whose presence may produce 
a significant perturbation on read-out fidelity. Thus the 
current pulse method is also sensitive to large splittings 
along multiple energy branches. 

MRT degeneracies also have very narrow bias value 
widths. For example, the inset of Fig. 2(a) shows a split- 
ting width of less than 0.1 nA. The bias widths become 
only smaller for splittings at larger Nl- The horizontal 
axis of Fig. 2(a) corresponds to more than 300 nA. Typ- 
ical experiments sampling only a limited number of bias 
values likely probe only a subset of the (many) MRT split- 
tings. Changes in experimental conditions (e.g. bias drift 
and noise, or temperature cycling) may generate surpris- 
ingly large shifts in the observed splitting distributions 
if they result in a different subset of sampled MRT de- 
generacies. These or other features could result in transi- 
tion spectra with a varying distribution of splitting sizes 
and bias-value locations which, due to their complexity 
and variability in time, might appear to have a micro- 
scopic origin. Such variations seem more consistent with 
a model of microscopic critical current fluctuators, sug- 
gesting that both MRT and microresonator effects are 
present If this is the case, it is important to identify 
which observed splittings are due to which mechanism. 

Measured with sufficient resolution, the transition fre- 
quency avoided crossings due to MRT should have dis- 
tinctive characteristics. When driving — > 1 transitions, 
a splitting in the |0) £ branch should produce crossings 
like that shown in Fig. 3(a), whereas a splitting in the 
|1) L branch should produce crossings like that shown 
in Fig. 3(b). The observed shapes may be strongly de- 
pendent upon the experimental measurement technique. 
Bias noise could smear out the splittings in the horizon- 
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FIG. 3: (a) The distinctive shapes of avoided crossings in 
the measured transition frequencies for splittings in the lower 
branch, (b) The avoided crossing transition shape for split- 
tings in the upper branch, (c) The figure shows that Ar ~ Al 
over an extremely large range of double-well circuit parame- 
ters. Bold lines show the splitting magnitudes A along the 
til = 0, 1, and 2 left well energy branches, with f3 = 4.5, 
N L = 3, and I c = 10/^A. 



tal direction. For splittings in the lower energy branch 
[Fig. 3(a)] this would leave a distinct frequency gap, 
givings observed splittings horizontally smeared appear- 
ances like those observed in Q, Q. In contrast, it is 
unclear if splittings in the upper branch [Fig. 3(b)] are 
consistent with observation. Improved experimental res- 
olution that revealed these distinctive avoided-crossing 
shapes would be compelling evidence for MRT. Han et al. 
have explored other complexities that arise when measur- 
ing systems that exhibit MRT 0. 

We derive an analytic expression for the energy split- 
ting between pair-wise degenerate left and right states in 
an asymmetric double well in the WKB approximation 
[l3|. This yields the splitting formula 
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7T n L \m R \e n ^+ m R+ 1 



(5) 

where S = J£ ^/2m [E - V ( 7 )]d 7 , m = C7($ /2tt) 2 , 
7i,2 are the classical turning points for the barrier given 
by ^(71,2) = E riL , and A L ~ huj L , A R ~ hu R are the 
left and right well level spacings at energy E nL . For deep 
right wells, Eq. iJHJ becomes independent of m R . In this 
limit, together with the cubic approximation accurate for 
shallow left wells, the splittings are approximately 
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(6) 



For the right well level spacing we use the WKB estimate 
A^? = 2wh/T c i 14], where T cl is the classical period of 
right-well oscillations with energy E nL . Splittings calcu- 
lated from Eq. © are shown as solid lines in Fig. 2(d). 
The agreement with the exact splittings (solid points) is 
excellent for lower lying states, and surprisingly good for 
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the excited states. Note that the tunnel splitting for- 
mula in Eq. predicts splittings exponentially larger 
than continuum tunneling rates: A sp ii t ting/rtunneiing ~ 
exp (18Nl/5) , making MKT effects important even when 
continuum tunneling is negligible. 

We have compared MRT splittings with Eq. for a 
number of numerical examples with Nr ~ 100 — 1000, 
but in principle one can fabricate circuits with many 
thousands of right- well states. The WKB formula for 
the splittings and level spacings allows analysis of cir- 
cuit parameters for very deep right wells where numerical 
treatment is impractical. Figure 3(c) shows Hlvl — Al 
(dashed line) and Ar = 2irH/T c i (thin-solid line) versus 
the ratio E c /Ej for I c = 10/xA, N L = 3, and = 4.5 
just below the threshold where the potential develops 
three wells. (For the circuit parameters in Fig. 2 and 
2], E c /Ej ~ 10~ 4 - 10~ 6 .) The value of I c determines 
the frequency scale on the left of Fig. 3(c) but leaves 
the relative positions of the plotted lines essentially un- 
changed. The top axis shows Nr from the harmonic os- 
cillator approximation. Observe that, perhaps unexpect- 
edly, Ar ss Al even for extremely asymmetric double 
wells. The bold-solid lines show the WKB splitting A 
when til = 0, 1, and 2. The validity condition for MRT 
A -C Arx is satisfied over a large range of circuit pa- 
rameters, and for Nr ~ 10 5 and greater. 

Dissipation suppresses resonant tunneling when Tr > 
An, where Tr ~ NrH/Ti is the width of excited right 
well states, and T% is the dissipation time for \\)r — > 
|0) fl 0, E|. Using the WKB expression for A R , we 
find the condition Nr < lu^Ti for observing MRT. For 
a phase qubit with w l /2tt ~ 10 GHz and T x ~ 10 - 100 
ns, resonant tunneling should be detectable as long as 
Nr < 600 — 6000 states. For the circuit parameters in 
Fig. 2 Nr ~ 600-700 and for those in 2] N R - 150-300, 
with a measured Ti ~ 25 ns. Thus, we do not believe that 
dissipation will remove the effects of MRT. If the intrinsic 
dissipation is actually much smaller so that T < A id, 
it should be possible to observe coherent oscillations ljj . 

In conclusion, we show that significant MRT effects 
should be present for extremely asymmetric double well 
phase qubits, and thus MRT should be taken into ac- 
count in the important effort to fully characterize mi- 
croresonators or other splittings mechanisms. Our anal- 
ysis provides tools and can guide experiments to help dis- 
tinguish between three main possibilities: (1) Both MRT 
and microresonators are present, (2) MRT effects explain 
all the observational data, and (3) MRT is entirely ab- 
sent. We believe that (1) is most likely; however, due to 
the complexity of effects from MRT, further experiments 
and detailed modeling are necessary to definitively rule 
out (2) and (3). Finally, our Letter provides general tools 
for exploring the quantum mechanics of extremely asym- 
metric double- well systems. 
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